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We describe a variational calculation for the problem of screening of a point charge in a layered 
correlated metal for dopings close to the Mott transition where the screening is non-linear due to 
the proximity to the incompressible insulating state. We find that external charge can induce locally 
incompressible regions and that the non-linear dependence of the screening on density can induce 
overscreening in the nearest nearby layers while preserving overall charge neutrality. 



INTRODUCTION 

Correlated metals can exhibit a large variation of elec- 
tronic compressibility for small changes in the electron 
density. An external charge can have varying conse- 
quences depending on this compressibility. This may 
lead to a large change in the local compressibility, to the 
point that the metal may be driven into a local insulat- 
ing, incompressible state near the external charge. Such 
non-linear screening is further complicated by the layered 
structure of the interesting correlated metals, such as the 
cuprates, since an additional length scale, the layer spac- 
ing is introduced. We demonstrate that these aspects of 
the problem can be analyzed quantitatively by a vari- 
ational calculation which gives further insight into the 
physics of the problem. The current interest into the 
problem of electrostatic screening in a layered correlated 
has been induced by recent muon spin relaxation exper- 
iments in underdoped metallic cuprates. [1-8] 



SCREENING IN A LAYERED METAL. 

Electrostatic screening in a layered metal can be ana- 
lyzed by a model of an equidistant stack of metallic planes 
separated by a dielectric medium. We consider the prob- 
lem of a fixed charge in between the planes while the 
screening charge is on the metallic planes. We start with 
the Poisson equation, written in terms of electrostatic po- 
tentials, cf>(r,z n ), and screening charge, p(r,z n ), on the 
metallic planes, 

n' 

where r is a continuous planar position, z n = (1/2 + n)c 
is z-axis position of the metallic plane in the stack, c 
is distance between the planes. [10] Here M (r,z„) = 
(|e|/e)/-y/r 2 + (z n — d) 2 is the electrostatic potential of 
the external charge. The external charge is between 
the two planes, at z n = ±c/2, at position d above the 
center. Going to planar momentum variable, via 



<K£) = (1/2 71 ") J d 2 r exp{— i£r}(f>(r) we write 
</>{£, Zn) ~ <t>iu,(£,Zn) = 47r|e| 2^p(£, z n ) — (1) 
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Due to rotational symmetry the planar momentum £ en- 
ters only via £ = |£|. Here <^(£, z n ) = 4 7 r|e|e-«l 2 "- c V 2 £ 
is the potential of an external charge. Rewrite it in the 
form 

z n - 4) " z 'n)) = 47rep(£, z n ) 
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(2) 

where A = A~~ l is Laplacian operator on a finite stack 
of metallic planes defined as a matrix inverse of A nn i — 
e -i\z-n~z ,,| j2^_ For a finite stack of N planes it is a tri- 
diagonal matrix of the form, 
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where g' = £(1 + cothc£),g = 2£cothc£ and / = 
— £/ sinhc£. Eq. (2) can be obtained as a condition for 
stationarity of the functional 

nn' 

+ S F ({<Kr ) z n )}) (4) 

with (f> as an independent variable. Here z n ) — 
4>(£.i z n) — 4>u(£,i z n)- The screening charge density 
is defined via the second term by Air\e\p{r 1 z n ) — 
—dSp/dip(r,z n ). For a metallic plane characterized by 
linear compressibility we have p = — («/27r|e|)0 where 
k = (2ir \e\ 2 / e)(dn/ dp) is the inverse screening length of 
the two-dimensional metal [9] and 

S °f = j 'I^E 2 ^™) 2 (5) 
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The form of Sf for non-linear screening will be discussed 
later. We will take the simple view that for small dis- 
placed charge the metallic plane has a large and fixed 
charge compressibility and therefore the screening is lin- 
ear. Thus in the situation where most of the screen- 
ing of external charge is done by the two nearest planes, 
z n = ±c/2, the screening charge in the rest of metallic 
planes is small and therefore can be accounted within 
linear screening regime. In this situation one has to con- 
sider the effects of non-linear screening in the two nearest 
metallic planes only. 

To obtain a functional that depends on the potentials 
in the two nearest planes only [11] we rewrite Eq. (4) in 
the matrix form (the ^-integrals are implied) 
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We have kept the (linear) compressibility term for all 
planes except the two in the center, z n = ±c/2. Here 
\& = {ip{£, z n )} n is a vector planar index as an index. 
Similarly, $ M = {<j>n(^,z n )} n . We separate out the two 
nearest planes $o — {^(£> 2+1/2), Z-1/2)}) the rest 
of planes above, *ff + , and the rest of planes below, 

(and similarly for $„). The 2x2 block, A o, multiplies 

(2) 

the two planes near the external charge. The block A^g 
is 2 x ^y^- block with only one non-zero element in cor- 
respondence with the tri-diagonal structure of matrix A. 
Omitting v&o-independent terms, we write 



S = ^t[A 00 + Ai]* - + H*o] 

Ai = -A 0+ [2k + A ++ ]- 1 A +0 + (+++-) 

$ M = 2kA 0+ [2k + A ++ ] _1 $^ + (+++-) (7) 



which defines the bilinear and linear terms in the func- 
tional for two central planes. The first term in the expres- 
sion for Ai has a form [qq] with B = (— £/ sinhc£;) 2 [2ft+ 
A( ++ )]~ 1 | A r i 7v, (here N is the dimension of the A ++ ma- 



trix). We define matrix D via 2k + A 



(++) 



[2£ coth c£ 



2k]D. It is equal to 1 on the main diagonal and — l/(2a) 
on the two near diagonals, a(£) = coshc£ + sinhc£. 
Taking the limit N — > 00, ignoring a', and using the 
identity [£) _1 ]jviv = (1/2tt) d92a sin 2 9/[a-cos0] we 

obtain B = (£/ sinhc£) [a — y/a 2 — 1]. Adding the ( ) 

block in Eq. (7), 
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The linear terms in Eq. (7) are evaluated similarly, 
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We obtain (here <j>± 
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Yp~\M z n = Mk+y^)) (9) 

(z„ = c/2) ± 0(z„ = -c/2)), 

1 E{^± G ±^±+ 2< ^± F ±}+ / d2rS F 

^2 cosh c£ — a + \f 1 
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sinh c£ 
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sinh c£ 



sinh ^ cosh d£ 
- cosh ^ sinh d£ 



(10) 



This functional serves as a starting point for an analysis 
of the effects of non-linear screening in the two nearest 
planes assuming that all other metallic planes are in the 
linear screening regime. When the screening in the two 
nearest planes is linear, Sf = j[2k((/> + + <f>- )], Eq. (10) 
is solved with 



87r|e| 



sinh s cosh d£ 
- cosh y sinh ^£ 



+ "v/c^ 2 
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which can also be obtained directly from Eq. (1). The 
planar screening length is controlled by the singularity of 
(f>° (£ ) in the complex plane of £ which is closest the real 
axis, i.e., one of the branch points £* which is a solution 
of a(£*) = ±1. The screening length approaches c/n for 
kc > 1. 



NON-LINEAR SCREENING. 

To discuss the non-linear screening, we need to specify 
the compressibility part of the functional, Sf, and vari- 
ational functions. We take a simple model in which the 
non-linear compressibility in the metallic plane is a step 
function: when local screening charge density exceeds a 
certain threshold, p* , the metallic layer goes over into 
a local insulating (incompressible) state, 5p/6(f) — 0. In 
the cuprates, for instance, we have p* ~ 0.1 holes per 
unit cell around optimal doping. In the absence of ex- 
ternal charge, the local hole density is determined by 
chemical doping, p = 0. At high enough doping the 
copper-oxide plane is metallic and it is characterized 
by linear electrostatic response, 5p/8cf) = — /t/(27r|e|/e). 
We note that when the threshold density p* is small a 
significant fraction of the electric field may leak to the 
next plane, z n = ±3c/2, and be strong enough for it to 
reach the threshold density as well, which will invalidate 
our assumption of all other planes being in the linear 
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screening regime. The extension of Eq. (10) to allow 
non-linear screening in more than two nearest metallic 
planes is straightforward. Since the local charge den- 
sity is defined via Sp, this model requires Sp = K(f> 2 
for 4> < 4>* and Sf = K(/>*(2(f) - 4>*) for <j) > 4>* . Here 
2K<fi* = 47r |e |p* . The effect of non- linear compressibil- 
ity is that the screening charge is pushed out of the im- 
mediate vicinity of the external charge (when compared 
with the screening in the linear case), i.e., the screen- 
ing cloud is larger. It is important that this does not 
affect the screening cloud at large distances where dis- 
placed charge is small and the screening is linear. This 
observation leads to the choice of our variational func- 
tion: the electrostatic potential for a screening of a point 
charge in a non-linear metal is modelled with a screen- 
ing potential in a linear metal due to a "smeared" , finite 
size, external charge. Such variational function can be 
constructed starting with Eq. (11). We further assume 
that the charge is only distributed in the direction par- 
allel to the plane with charge density v{r), i.e., the vari- 
ational function is </)(£, z n —±c/2) — v(£)<f>o(£, z n =±c/2) 
where (j)o is given by Eq. (11). The simplest choice is 
a constant charge density over a disk of radius R, i.e., 
v(r) = Q9(R — r)/(nR 2 ) where 9{x) is a step function. 
For the situation discussed below a non-uniform distri- 
bution v(r) = Qcxpj— r / R} /(2irR 2 ) proves to be better. 
For the non-linear screening the density is not simply re- 
lated (not proportional) to the local screening potential, 
thus Q is an independent variational parameter. The in- 
terpolating property of our choice of variational function 
is that in the linear screening case it is an exact screen- 
ing potential with R = and Q = 1. We now discuss 
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FIG. 1. Non-linear screening of the point charge +1 in a 
layered metal. We take kc = 10 and c = 1.5a. The gray line 
is the non-linear screening density profile of an external charge 
in the center, d = 0. The blue and red squares represent the 
screening density in the two planes for asymmetric position 
of an external charge, d = c/4. The solid red/blue lines show 
the linear screening density profile for the same geometry in 
the upper and lower metallic planes. 

screening of an external positive charge, +|e|, in a lay- 



ered metal which geometry which roughly corresponds to 
LSCO cuprates, presented in Fig. 1. The planar distance 
is in units of planar unit cell size, a, and the density is 
given in units of charge per planar unit cell. For a sym- 
metric position of the external charge, the insulating re- 
gion has a radius of 0.68a, and it contains about —0.30 
charge (in units of +|e|). Note that the for large val- 
ues of linear compressibility (kc > 1) the density outside 
the insulating region decays rather slowly, at distances of 
about inter-plane distance, c, and the charge density de- 
crease to a half of its value in insulating region at about 
la from the center. The total charge in each plane is 
—0.5, ( —0.48 for linear screening in the same geometry). 
These numbers support our assumption that all other 
metallic planes are in the linear screening regime. The 
situation changes slightly if we shift the external charge 
away from the symmetric position, d = c/4, which more 
closely reflects the muon site in the LSCO lattice. The 
area of the insulating region in both planes increases in 
this asymmetric situation, = 0.69a in one plane and 
1.05a in the other. The total screening charge in each 
plane is, —0.58 and —0.74 respectively (—0.26 and —0.69 
for linear screening in the same geometry). Compared 
to the symmetric external charge position, in this asym- 
metric situation the plane which is closer to the external 
charge is in the locally insulating state over a larger area 
and therefore, the effective, partially screened charge seen 
by the other, more distant plane, is also larger. The to- 
tal charge of the two planes, —1.33, is greater than one, 
i.e., in this geometry the non-linear screening leads to an 
overscreening of the external charge in the two nearest 
metallic planes. The integrated screening charge in the 
infinite stack planes remains equal to the external charge. 

We now compare our results for the non-linear screen- 
ing in LCSO to a recent dynamical mean field calculation 
presented in Ref. [8]. We consider the non- linear screen- 
ing for a symmetric position of the external charge. Our 
variational method gives a change in the local density of 
about 0.10 holes per (planar) unit area in the first cell. 
Ref. [8] gives a change of about 0.07 holes per planar unit 
area. In the adjacent cells, our variational method and 
Ref. [8] give a change of about 0.01 holes per unit area. 
Both models give negligible changes in the local density 
of holes for cells located beyond the adjacent cells. 

CONCLUSIONS 

We have developed a simple variational framework for 
the analysis of the screening of point charge in the cor- 
related layered metals, such as cuprates in the metallic 
doping range. Applied to the screening of the muon in 
the LSCO cuprates it supports an earlier qualitative dis- 
cussion in Refs. [7 and 8]. 

We point out that the method presented here may be 
used to study the charge distribution, rectification, tran- 
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sistor action and other surface effects between weakly 
doped Mott insulators and metals. 
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